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Abstract 

We study the vibrations of a hanging thin flexible rod, in which the dominant restoring 
force in most of the domain is tension due to the weight of the rod, while bending elas- 
ticity plays a small but non-negligible role. We consider a linearized description, which we 
may reduce to an eigenvalue problem. We solve the resulting singularly perturbed prob- 
lem asymptotically up to the first modification of the eigenvalue. On the way, we illustrate 
several important problem-solving techniques: modeling, nondimensionalization, scaling, and 
especially use of asymptotic series. 



1 Introduction 

Linear transverse vibrations of a uniformly stretched string are modeled with the wave equation 

PWtt - TWyy = 0, (1) 

where w(y,t) is the displacement of the string, p the mass per unit length, and T the tension. In 
([lj the string is assumed perfectly flexible; if resistance to bending is not negligible (e.g., for a 
rod), then the equation requires an additional, fourth-order, term, 

PWtt - TWyy + ElWyyyy = 0, (2) 

where E is Young's modulus and / is the cross-sectional moment of area. If the tension depends on 
position y, then the middle term of ^ is replaced by a variable-coefficient operator in divergence 
form; in particular, if the rod is vertical so that the tension at a point y results solely from the 
weight of the rod below it, located say in [0, y], then the tension is T(y) — pgy and ^ is modified 
to 

PWtt - Pg(yWy)y + ElWyyyy = 0, (3) 

where g is the acceleration of gravity. In this paper we study ^ in the asymptotic limit where 
the bending forces are small compared to tension, using singular perturbation^] 
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Our interest in this problem originated from Manela and Howe's paper (MH09| simulating the 
waving of a flag in the wind. The displacement of the flag was represented as an expansion in 
terms of eigenf unctions of the second-order linear operator 



d ( dw 
dy \ dy 



) 



i.e., the part of |3]) due to tension alone, neglecting bending resistance. It was found that this 
expansion converged very poorly, if at all. In this paper we study the effect of a small bending 
resistance on the eigenf unctions, without reference to wind-driven motion. 

The problem of fluttering plates and rods driven by a fluid has been studied earlier in various 
places. In a recent paper, Argentina and Mahadevan [AM05 study the problem where the bending 
rigidity dominates, and in a previous paper, Dowling Dow87 uses matched asymptotic expansions 
to solve a slightly different problem where the fluid- loading is taken into account. This has the 
consequence of moving the main singular point into the bulk which results in significantly different 
behavior. 

2 Preliminary Analysis 
2.1 Boundary Conditions 

Equation ([3| must be supplemented with two boundary conditions at each end of the domain, say 
< y < L. At the free end y = there should be no bending moment (i.e., Elw yy (0,t) = 0) and 
no force. At a point where the tension vanishes, the force is given by the negative derivative of 
the bending moment, so we obtain the second boundary condition w yyy (0,t) = 0. 

At y = L we consider two distinct possibilities for boundary conditions: clamped and pinned. 
Both imply that there is no deflection at y = L; thus, w(L, t) = 0. For the second boundary 
condition, clamped imply zero slope (i.e., w y (L,t) = 0), and pinned imply zero bending moment 
(i.e., w yy (L,t) = 0). 

Thus, in summary, the two possible sets of boundary conditions are 



2.2 Non-dimensionalization and Scaling 

Despite there being several parameters that control the behavior of pi), these may be reduced to 



a single non-dimensional "group." Let us scale y by the length L and t by L/c, where c = y/gL 

specifies the order of magnitude of the speed tension waves; i.e., let 



The reason we use the timescale of the tension waves is that we are interested in looking at the 
thin rod case, where the elasticity plays a small role in the dynamics. (We do not scale w since in 
a linear equation this would not change anything.) 




w y (L,t) = 0, (Clamped) 
w yy (L,t) = 0. (Pinned) 



(4) 





(5) 



Substituting (|5j into ^ yields 



w tt ~ {y w v)y + £w yyyy = °> 



where e = 



EI 



(6) 



pgL 3 ' 
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The dimensionless constant e compares the importance of bending elasticity (0(EI/L 2 )) with the 
maximum tension (pgL). Of course e is small if L is large. Let us also examine the dependence of 
e on a, the width of the rod. The second moment of area is given by 

x 2 dxdz, (7) 



c 

where C is the cross-sectional area of the rod and x is the direction of bending. This moment 
scales like a , while p, mass per unit length, scales like a 2 . Thus, e contains an implicit factor of 

0=0 

Below we omit the tildes from ([6]). 
2.3 Separation of Variables 

We look for a solution of ^ in separated form w(y, t) = u(y) x f2(f), and find that u must satisfy 

eu"" - {yu')' =\u, < y < 1, (8) 

where A is the eigenvalue parameter, and 

n(t) = fl e ±lVJt . (9) 

Note that ^ is singular for two reasons: (i) e multiplies the highest-order derivative in the 
equation and (ii) the coefficient of leading-order derivative in the reduced equation (after setting 

£ = 0), 

- (yu 1 ) 1 = Xu, (10) 

vanishes at one end of the interval. In the remainder of the paper we solve asymptotically this 
singularly perturbed eigenvalue problem, subject to the boundary conditions 

f either = 0, (Clamped) 

u"(0) = 0, u"'(0)=0, «1) = 0, plus „ ' \ ' (11) 

I or u (1) = 0. (Pinned) 

Since eigenfunctions are determined only up to a multiplicative constant, we add the normalization 

u(0) = 1, (12) 

thereby selecting a unique solution of this problem. 

Incidentally, we claim that ^ with either boundary conditions (11) is self-adjoint and in fact 
positive-definite. This may be proved with the usual integration-by-parts argument, with one 
subtlety: usually, boundary terms vanish because of the boundary conditions, but the term (yu')' 
makes no contribution at y — only because the coefficient y vanishes there. Thus, all eigenvalues 
of this problem will be positive real. 



tin |3|, the coefficient EI of the fourth derivative is appropriate for a solid rod but not for a cable or string — i.e., 
many small fibers twisted together. The bending resistance of any one fiber in a string is all but infinitesimal; the 
primary resistance comes from the friction of fibers sliding over one another. While it is difficult to calculate the 
bending resistance of such a collection of fibers, this resistance is much smaller than the Young's modulus times the 
area moment of the whole cable. Thus, appropriate values of e for a string may be very small indeed. For accurate 
modeling of a string it might be necessary to include a friction term, say proportional to wt, in pf. Such a term 
would not change the eigenfunctions found below, and its effect on the time dependence |9jl is easily calculated. 
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The remainder of this paper is organized as follows. In Section 3] we propose a naive derivation 
of the limiting behavior of the eigenvalues and eigenfunctions of (18 ) as e — > 0, and we present some 
numerical results. In Section [4] we find the boundary-layer scalings for the asymptotic analysis 
and solve the reduced equations in various regimes. In Sections [5] and [6] we derive matched asymp- 
totic series for the eigenvalues and eigenfunctions with clamped and pinned boundary conditions, 
respectively. Finally, in Section [7] we discuss the applicability of small-e asymptotics for later 
eigenvalues in the sequence of eigenvalues. 

In addition to its research interest, this problem provides a relatively simple example of a 
matched asymptotic expansion that requires logarithmic terms. Thus for pedagogical reasons, we 
strive for careful, thorough explanations, and we have included several exercises for the readers to 
sharpen their understanding. 



3 First Attempts 



3.1 The Naive Solution: e = 



Setting £ equal to zero in ^ yields equation ( 10 1. Of course no solution of this equation can satisfy 
all the boundary conditions (11). It is natural to conjecture that at y = 1 the "most physical" 



solution will at least satisfy the lower-order boundary condition there 

= 0. 



(13) 



Because (10 1 is singular at y = 0, it is unclear what boundary conditions, if any, ought to be 



imposed there. 



In fact, (10 1 can be solved explicitly in terms of Bessel functions of order zero. The connection 



to Bessel functions may be motivated by expanding all terms in (10), 

yu" + u' + Xu = (14) 
and observing that this equation, like Bessel's equation, has exactly two singular points, a regular 



one at y — and an irregular one at infinity. Moreover, the indicial equation of (14), obtained by 
seeking a series solution 



(15) 



where Co ^ 0, has a double root p = 0, again like Bessel's equation of order zero. To make the 
reduction, consider the substitution u(y) = v(Cy p ); a simple calculation shows that if one chooses 
p = i C = 2\J A, then v satisfies Bessel's equation v" + x~ 1 v' + v = 0. Thus, the general solution 



of ( 10 ) is 



u{y) = aJ (2v/Ay) + bY (2^) 



(16) 



for arbitrary constants a and b. As y — > 0, Jq is smooth, but Y blows up logarithmically (see 



Appendix A.2). Unless a = b = 0, neither of the boundary conditions (|TT|) at y — can be 



satisfied. Let us require that b — so that at least u remains bounded (and \12\ is meaningful). 
If the boundary condition ( 13 ) at y = 1 is to be satisfied, then A must satisfy 




Ji 



o (2VA) 



0. 



(17) 



This intuitive discussion gives eigenvalues A w 1.4458,7.6178,18.721,... with eigenfunctions 
Jo(2y/Xy), for either clamped or pinned boundary conditions. The eigenfunctions are graphed 
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in Figure [T] Despite the many loose threads in the argument, exactly these eigenvalues and 
eigenfunctions will emerge as the leading-order term in our asymptotic solution, thus providing a 
far more satisfactory derivation. Moreover, the difference between clamped and pinned BC will 
emerge in the higher-order terms of the series. However, before tackling the asymptotics, we turn 
to numerics. 




A = 1.4458 
A = 7.6178 
A = 18.721 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

V 



Figure 1: The naive solution of the reduced equation (10 1 using the three smallest roots of the Bessel 
function Jo(2vA). 



3.2 Numerical Solutions 

In an exploratory numerical code, we approximated (|8| by differences on a uniform grid. However, 
especially for small s, we found this simple approach gave unreliable results, even when up to 50,000 
grid points were used. For the data presented below, we used the boundary- value solver bvp5c 
in Matlab. The interval < y < 1 was divided into 3 distinct subintervals (corresponding to the 
2 boundary layers and the bulk domain, introduced below), which were connected with internal 
boundary conditions. Results from the uniform-grid code were used as initial guesses for the 
internal iterative solver. Very small values of e were approached by continuation. 

Graphs of the computed eigenfunctions resembled the naive eigenfunctions in Figure [T] but 
graphs of their derivatives differed substantially. This issue is explored in some detail in Subsec- 
tions [5^4] and [fT^e) ■ 

Figure [2] shows the divergence of computed eigenvalues from the naive approximation; specifi- 
cally, a log-log plot of |A(")(e) - AW(0)| /A(")(0) vs. e, where AW(0) is the n th root of J a (2V\). 
The results support the naive analysis and also suggest that 

A(")(e) = A (n) (0) +0(e p ) (18) 

where p — 1/2 or p = 1 for clamped or pinned boundary conditions, respectively. This suggestion 
is confirmed by the asymptotic series solution constructed below. 
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(a) Clamped 



Figure 2: A log-log plot of the relative error e(e) 



A(">(e)-Al"'(0) 
A<")(0) 




logio e 

(b) Pinned 



Mi 



for the first five eigenvalues. Accord- 
ing to our asymptotic analysis, in the clamped case, e(e) ~ (A 1 / 2 /Ao)e 1 ^ 2 where \i/ 2 is given in Table^ 
since \\/2 = Ao, all lines collapse onto e(e) = e 1,/2 (which is drawn for reference). In the pinned case 
e(e) ~ (Ai/Ao)e where Ai is given by (111 I. Reference lines e(e) = (Ai/Ao)e are drawn in the figure. 



4 Preparation for the Asymptotic Solutions 
4.1 Boundary-Layer Scalings 

In the bulk of the domain, far away from the boundaries, we assume that e times the fourth 



derivative is small enough that we can ignore this term to lowest order, obtaining (10). As 
we have seen, solutions of this equation cannot satisfy all the boundary conditions. We expect 
therefore that there are two adjustment zones — boundary layers — connecting the solution in the 
bulk to the boundaries, in a way that satisfies the boundary conditions. Near the boundaries, we 
expect that the fourth derivative becomes large enough that, even when multiplied by e, it cannot 
be ignored. As a first step we have to find the boundary-layer scalings that are appropriate near 
the boundaries. 

To find the scaling near the y = boundary, we assume that y is of order e p ; specifically that 
u(y;e) can be approximated by U(X;e) where X = s p y. To find p wc differentiate according to 
Eq. (|8| and find that 

s 1+4p U"" - e p (XU')' = XU (19) 

a be 

where prime indicates differentiation with respect to X. An appropriate scaling must balance 



Eq. (19); i.e., two terms must have the same order in e and the third term must have order this 



high or higher. Thus, we look at three possible cases: 

a«i Thus 1 + 4p = p and so p = —1/3; in this case terms (a) and (b) are much larger (i.e., lower 



G 



order) than (c). Hence it is a good scaling of the boundary layer. 

b « c Thus p = (i.e., no scaling happens.) Here the two terms (b) and (c) are larger than (a), 
so this scaling is valid; indeed this is the bulk scaling. 

a sa c Thus l+4p = 0sop — ~ 1/4 and so term (b) is much larger than (a) and (c), so this scaling 
does not balance]!] 

The only boundary layer we find near y — scales the inner variable X like e 1 / 3 , and the resulting 
equation is therefore 

U"" - {XU'Y = s 1/3 XU. (20) 

To find the scaling near at the y = 1 boundary, we write the solution u(y\ e) as V(Z; e) where 
Z = e q (l - y). From Eq. ^ we obtain 

£ l+ig v m> _ £ 2 q ^ _ £ - q ^ y,y = w (21) 
Here there are also three possible cases. 

Exercise 1. In the same way as before, show that the only valid, interesting scaling for small e 
at the y = 1 boundary has q — —1/2, and that the resulting equation for V(Z; e) is 

V"" - ( (l - e 1 ^) 1/') ' = eXV. (22) 
4.2 Reduced Equations 

In the asymptotic series for the solution, terms of order 1/3 and 1/2 are forced by the boundary 
layers, and subsequent terms include all sums of integer multiples of 1/3 and 1/2. Thus, we look 
for a solution that has the following form: 

• In the bulk we expand the solution u(y; e) as 

u(y; e) = u (y) + e 1/3 u 1/3 (y) + e 1/2 u 1/2 (y) + e 2/3 u 2/3 (y) + e 5/6 u 5/6 (y) + e Ul (y) + ... (23) 

• Near the y = boundary, the solution is approximated by a similar asymptotic series for 
U(X;e), with X = e^^y. 

• Similarly, for the solution near y = 1 with V(Z;e), where Z = e~ 1//2 (l — y). 

• The eigenvalu^]\ is itself expanded as an asymptotic series: 

X(e) = X Q + e 1/3 A 1/3 + e 1/2 A 1/2 + e 2/3 A 2/3 + e 5/6 A 5/6 +eX 1 + ...s (24) 

Below, we will see that the bulk series needs to be augmented by a term of order e log e, and 
logarithmic terms will be needed at higher order as well. However, for now we wait for this 
complication to arise naturally. 

tThis scaling would be valid in the opposite asymptotic limit where one considers e to be large. 
§Of course, there is a sequence of eigenvalues with corresponding eigenfunctions. We drop the superscript in the 
\( n ) notation for sake of a cleaner presentation. 
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For the three domains, we substitute the asymptotic expansion of u, U , or V, and of A into 
the appropriate differential equation and collect powers of e. Doing this for an expansion up to e 1 
results in 



for u, 



e° 




X uq = 






£ l/3 


{v u 'i/ 3 y 


+ Ao"l/3 


-A1/3U0 




£ l/2 


(yu' 1/2 )' 


■ A Mi/2 


-A1/2U0 




£ 2/3 




' A M 2 /3 = 


-A1/3U1/3 - 


- A2/3M0 


£ 5/6 


(y u 5/ e y 


+ A0M5/6 = 


-A1/3U1/2 - 


- A 1 / 2 u 1 / 3 - 


e 1 


(y<)' + 


X u 1 = Uq" 


— A 1 / 3 U 2 /3 


- A1/2U1/2 



(25) 



A2/3 M l/3 - AxU 



for U, and 





77"" 


(XUtf = 







£ l/3 


77"" 

U l/3 




= A E/ 




£ l/2 


77"" 




= 




£ 2/3 


77"" 

U 2/3 


- (xu^y 


= MUl/3 


+ A1/3C/0 


£ 5/6 




(xu' 5/ J 


= A0C/1/2 


+ Ai/2^0 


e 1 


U{"' 


(my = 


A0U2/3 + 


Al/3^1/3 + A 2 /3^0 



(26) 





V"" 


- V" - 


-- 


£ l/3 


V"" 


— V" 
V l/3 




£ l/2 


V"" 
V l/2 


— V" 




£ 2/3 


V"" 

V 2/3 


— V" 
V 2/3 




£ 5/6 


V"" 
V 5/6 


- n' /6 




e 1 


V"" 


- V{' = 


A 



(27) 



<zv; /3 y 

Vo - (ZV( /2 ) 



for V. 



4.3 Solutions of the Reduced Bulk Equations 



In the bulk, the reduced equations (25 1 have the form (yu')' +\qu — g, or expanding the derivative 
and dividing by y, 

„ 1 , A g{y) 



-u H u = 

y y y 



(28) 



We saw in Subsection 3.1 that in the homogeneous case, g = 0, the general solution of ([28]) is 
given by (16). To shorten the notation and form a more convenient basis for the solution space, 
we define 

J(y) = J {2^y) (29) 

y(y)=7ry (2 v ^o2/)-(logA + 2 7 ) J(y) (30) 

where 7 is Euler's constant. The modification of Yq simplifies the asymptotic behavior as y — > 
(see Appendix A. 2): 

Y{y)=\og{y) + 0{y\ogy). (31) 

Note that this notation hides Ao, an unknown constant that still needs to be determined. 
In the inhomogeneous case we have 
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Lemma 1. Given a continuous function g(y), there exists exactly one solution of (28), say w(y), 
that is continuous on the closed interval [0, 1] and satisfies w(0) = 0. 

For this singular problem, just one boundary condition suffices to determine the solution 
uniquely. 



Proof of Lemma^ (Existence) Using variation of coefficients, we find that one solution of (28 1 
is given by 



w(y) = ci(y)J(y) + c 2 (y)Y(y), 



(32) 



where 



ci(y) 
C2(y) 



1 



-Y(x)g{x)dx 



y i 

J(x)g{x)dx 



Y(x)g{x)dx 



o xW(x) 



J{x)g(x)dx. 



(33) 
(34) 



Here W = JY' — YJ' denotes the Wronskian of J(y) and Y(y), and we have used the fact 
W(y) — 1/y (see Appendix A.2). Regarding the first term in (32), although the integrand defining 
Ci in (33) diverges logarithmically at x = 0, the integral converges for all y £ [0, 1] and vanishes 
if y = 0. Regarding the second term (34), although Y(y) blows up logarithmically as y — > 0, the 



coefficient satisfies c^iy) — 0(y) so the product C2Y vanishes at y = 0. Thus, (32) solves (28) and 
satisfies the boundary conditions in Lem ma [Tj 

(Uniqueness) The general solution of (|28| has the form 

u (y) = w (y) + a J(y) + bY (y) 

for arbitrary constants a and b. For this function to be continuous we need b — 0, and for it to 
vanish at the origin, we need a = 0. □ 



Corollary 2. // Ao satisfies J Q (2\/X^) = and g(y) — J(y), then the solution to Eq. (28), w(y), 
does not vanish at y — 1. 



Proof. Since J(l) = 0, we have from (32) that 



;(1) = c 2 (l)Y(l) = / J 2 (x) dx ? 0. 

Jo 



(35) 



Of course the integral is positive, and Y cannot vanish at y = 1 since J already vanishes there 
and the pair's Wronskian is non-zero. □ 



4.4 Solution of the Reduced Boundary Layer Equations 



Equations (27) for the boundary layer at y = 1, which have constant coefficients, do not require 
any special discussion. Therefore, we focus here on the boundary layer at y = 0. 



The homogeneous versions of equations ( 26 ) have the form 



U" 



{XV 



(36) 
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We need to find four linearly independent solutions of this equation. By inspection, = 1 is 

ne W = U' , we n 

[W" - XW] = 0. 



one such solution. Observe that, if we define W = U' , we may rewrite (36) as 

d 



dx 



(37) 



Now W" — XW = is Airy's differential equation, whose solution space is spanned by the Airy 
functions, Ai and Bi. (See Appendix A.l for the definition and some elementary properties of 
these special functions.) Thus, 



U {2 \X) 



x 



Ai(x)dx, U (3) (X) 



x 



Hi(x)dx 



(38) 



provide two more linearly independent solutions. Since Bi grows super-exponentially as X — > oo, 
the solution cannot be matched to any solution in the bulk. By contrast, Ai decays super- 
exponentially so the integral to infinity converges; in fact, by (1201 



lim U (2) (x) = 1/3. 

X— Yea 



Incidentally, for use in the boundary conditions (11 121), we claim that 



t/ (2) (0) = 0, U {2) "(0) = Ai'(0) ^ 0, U {2) '"(0) = 0. 



(39) 



(40) 



The first relation is trivial; the second may be derived by differentiating (38) twice; and third may 

-1, 



be derived by differentiating (38) thrice and invoking Airy's differential equation. 



To obtain a fourth independent solution, we satisfy (37) by requiring that W" — XW ■ 
or since W = V , 

U'" — XU' = -1. (41) 
Specifically we let = ^ where ^ satisfies the following 
Lemma 3. There is a unique solution ^(A) of ( |41[ ) such that 

(a) *(0)=0, (b) ^"(0) = and 
(c) *(A) = logpO + 0(1) as X^oo. 

This lemma is proved in Appendix[B] and Figure [6] shows the graph of ^P, obtained numerically. 
Incidentally, the numerics indicate that as X — > oo 



*(A) = log(A) + * 00 +o(l) where 



1.3556, 



while it follows from (41 ) that 



*"'(0) 



-1. 



(42) 
(43) 



Below wc also solve inhomogeneous versions of equations (26). 



5 Asymptotics for the Clamped Case 

In this section, we calculate terms in the asymptotic series for the clamped case through the first 
non-trivial correction to the eigenvalue: i.e., order 1/2. The results are summarized in Table fl] 



10 



a 


A Q 


u a 


u a 


V a 





A 


J(y) 


1 





1 

3 








-A X 





1 
2 


A 


-X w(y) 





J'(l) [l-Z-e"*] 



Table 1: Coefficients of the asymptotic expansion for the clamped case. The leading-order eigenvalue Ao 



is a root of the equation Jo(2\/Ao) = 0. J and w are defined in (29 1 and (581, respectively. 

5.1 Zeroth-Order Solution 



(a) Boundary-layer near y — 0: As discussed in Section 4.4 the solution of the e°— equation 
in (|26|) is a linear combination 



U Q (X)=A + B f Ai(x)dx + CV(X), (44) 
Jo 

the integral of Bi having been excluded as unsuitable for matching. To satisfy the boundary- 



conditions (11 12), we require 

U (0) = 1, C / o(0) = 0, CC(0)=0. (45) 

By substituting ( |44| ) into the boundary conditions and recalling the derivatives of J Ai and ^ at 
X = 0, we deduce that A = 1,_B = C = 0, that is, 

Uo(X) = 1. (46) 
f'&j Boundary-layer near y = 1; The four-dimensional solution space of the e° equation in 



(27) is spanned by l,Z,e , e . Excluding e as unsuitable for matching, we write Vq(Z) 



A + BZ + Ce . The boundary conditions ( 11 ) for the clamped case 



Vb(0) = F '(0) = 0, (47) 
allow us to express two of the three arbitrary constants in terms of the third, yielding 

V Q (Z) = A(l-Z-e- z ). (48) 
The undetermined coefficient A will be found by matching with the bulk solution. 



(c) Bulk: We already know that the general solution of the £°— equation (25) in the bulk is 

u (y)=aJ(y) + bY(y). (49) 
The two constants will be determined in matching. 

(d) Matching: To match the bulk solution with that of the boundary layers, we compare the 
"outer limit" of the inner solutions (the boundary-layer solutions) to the "inner limit" of the outer 
solution (the solution in the bulk). Thus, near y = we need, ase->0, 



U 



(e^X) - Uo(X) = o(l) 



(50) 
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for an appropriate range of X. Specifically, we require that there exist numbers p, q, with < p < 
q < 1/3, such that (50) holds for all X such that 

e- p <A (51) 



In this case we may take p = 0, q = 1/3; i.e., the maximal range. It follows from (49) that 
u o(y) = a + blogy + 0(y\ogy) for small y. Hence, since e x l 3 X <§C 1, 



U 



(e 1/3 X) = a + b\og{e 1/3 X) + o(l). 



(52) 



Unless 6 = 0, the logarithm term in (52) is large and moreover depends on X. Therefore for uo to 
match onto Uq(X) = 1, we need a = 1 and 6 = 0. 
Near y = 1 we need 



u {l-£ 1/2 Z)-V {Z) =o(l) for Z in a range e~ p < Z < e~ q 



(53) 



where < p < q < 1/2. We take p = 0, q = 1/2. Given a, b as above, u (l - e 1/2 Z) = J(l) + o(l) 
while Vo(Z) = A(l — Z) + o(l). Thus matching the linear term in Vb requires that A = 0, and 
then matching the constant terms requires that J(l) = 0. Finally, extracting the implicit Ao from 
the argument of J, we obtain the characterization of Ao 

J {2^X~o) = 0. (54) 
This confirms the ad hoc solution we found in Subsection |3.1| and it verifies the first row of Tablc[T] 

5.2 The e 1/3 Correction 

Not much happens at this order, but it serves as practice for later calculations. 

(a) Boundary-layer near y = 0: At order e 1 / 3 , equation (26) is inhomogeneous with the 
right-hand-side XqUq = Xq, which has the particular solution U p (X) = — XqX. Ui/ 3 must satisfy 
boundary conditions analogous to (45 1, except that now Ui/ 3 (0) = replaces the condition Uo(0) = 
1. Since U p already satisfies the boundary conditions, we have that 



U 1/3 (X) = -X X. 



(55) 



As we shall see below, this term merely matches the first derivative of the bulk solution uq at 
y = 0. 



(b) Boundary-layer near y = 1; Arguing as in the e° case in Subsection 5.1 we deduce that 



V 1/3 (Z)=A(1-Z- 



(56) 



(c) Bulk: The e 1 / 3, equation (25) has the inhomogeneous term — Xi/ 3 uq. Thus, the general 
solution of this equation is 



ui/ 3 (y) = -Ai/ 3 w(y) + a J(y) + bY{y) 

where we define w as the solution (see Lemma [lj of 

(yw'Y + X n w = J such that w(0) = 0. 



(57) 
(58) 
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(d) Matching: Near j/ = Owe need, as e — > 0, 

(«o + £ 1/3 W i/ 3 )(e 1/3 X) - (U + e^ 3 U 1/3 )(X) = o(e^) 



(59) 



for X in a range e p <C X <C e 9 where < p < g < 1/3. We take p — 7 q — 1/6 so that 
(e 1 / 3 ^) 2 = o(e 1 / 3 ). Then by Taylor expansion at y = 0, 



uoie^X) = uo(0) + £ 1/3 ^(0)X + o(e 1 / 3 ). 



(60) 



Therefore, obtaining Up(0) from (126} in the Appendix, and recalling that w(0) = in ( p57} , we 
fincffl 

(«o + £ lA V/ 3 )(e 1/3 X) = 1 + e 1/3 (-AoX + a + blog{e 1/3 X)^ + o(e 1/3 ). (61) 
On the other hand, 

(U + e 1/3 C/ 1/3 )(X) = 1 - e x ' 3 \ X. (62) 
Thus (59) requires that a = b = 0. 



Near y = 1 we need 

(«o + £ 1/3 «i/ 3 )(l - e 1/2 ^) - (Vo + e 1/3 ^i/ 3 )(^) - o(eV3) 



(63) 



for Z in a range e p <C Z <C £ 9 where < p < q < 1/2. We take p = 0,q = 1/6 so that 
e 1//2 Z = o(e 1 / 3 ). Now by a Taylor expansion near y=l 



so by our choice of q 
Therefore, 



u (y) = -J'(i)(i~y) + O((i~y) 2 ), 
u {l-e 1 ' 2 Z) = O{e l ' 2 Z) = o{e 1 ' 3 ). 



(u + e 1/3 Ml/3 )(l - e 1/2 Z) = eV 3 u 1/a (l) + o(e i/3) = _ e V3 A + o(e V3) 



(64) 
(65) 

(66) 



where we have recalled that a = b = in (57). On the other hand, Vq = 0, so by (56 1, as Z — > oo 

(Vo + e 1/3 ^i/ 3 )(^) - As^il -Z) + oie 1 ' 3 ). (67) 

Thus matching requires that A = and A 1 / 3 u'(l) = 0. Recalling from Corollary [2] that w(l) ^ 0, 
we obtain A1/3 = 0, and we have verified the second line of Table [TJ 

5.3 The e 1/2 Correction 



(a) Boundary-layer near y = 0: Equation (26), which at order 1/2 is homogeneous, and the 
three homogeneous boundary conditions imply that t/1/2 = 0. 



(b) Boundary-layer near y = 1: Since Vq = 0, equation (27) at order 1/2 is homogeneous. By 
the same argument as for Vq and T4/ 3 , this equation and the boundary conditions yield 



V 1/2 (Z)=A(l-Z-e- z ). 



(68) 



^The term log(e 1 / 3 X) in this equation is a little disturbing. One expects the function that multiplies e 1 / 3 to 
depend on X alone, not e. This confusing behavior, a consequence of the singularity of ( | 1 [ > at y = 0, is not an 
issue here since the matching implies that 6 = 0. However, exactly this complication will force us to add a eloge 
term to the series in Section [6] below. 
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(c) Bulk: The e 1 / 2 equation of |25| ) has the same family of solutions as in order e 1 / 3 : 

ui/ 2 {y) = -X 1/2 w(y) + aJ(y) + bY(y) (69) 



where w is the solution of ( 58 ) . 



(d) Matching: Matching at y — as above, we deduce that a — b = in (69). 
At y — 1 we need 

(n + e^V/a + £ 1/2 u 1/2 )(l - e^Z) - (V + e^ 3 V 1/3 + e^ 2 V 1/2 )(Z) = o(e^ 2 ) (70) 

for Z in a range e~ p <C Z <C £~ g where < p < q < 1/2. We choose p = 0, q = 1/4 so that 
{e 1 ' 2 Z) 2 = o(e 1 / 2 ). Now by @ and g, 

(«o + £ 1/3 ui/ 3 + e 1/2 "i/ 2 )(l - e 1 ' 2 Z) = e 1 ' 2 [-J'{l)Z - \ 1/2 w(l)] + o^ 1 / 2 ), (71) 

while letting Z — > oo in ( [68] ) we deduce 

(Vb + e 1/3 F 1/3 + £ 1/2 ^ 1/2 )(Z) = e 1 / 2 A(l - Z) + oie 1 ' 2 ). (72) 



To match the expressions in (70) we need A = J'(l) and 

Ai/ 2 - -J'(l)Ml). 



(73) 



Remarkably, it follows from Lemma [4] that Ax/a = Ao- The verification of all entries in Table [T] is 
now complete. 

As a check on our calculations, in Figure [5b] below we present a log-log plot of the error in the 
two-term approximation A + A^e 1 / 2 to the clamped eigenvalues. The resulting lines of slopes 
near 1 suggest that the next non- vanishing correction to the eigenvalue will happen at the 0(e) 
level. An exercise at the end of Section [6] includes verifying this point. 



Lemma 4. The function w satisfies 



w{\) 



J'(l) 
Ac 



Proof. Recall formula (35) for w(l). Manipulating [AS64, 11.3.34], we find that 



J 2 {t)oIt= / J 2 (2^/\oT)dT = J 2 (2^/X ) + J 2 (2^/X ). 



(74) 



Since Ao is a root of Jo(2y/~Xo), we may drop the first term. Regarding the second we invoke [AS64, 
9.1.28] to obtain 

J 2 {2^Y Q ) = J l2 {2^X~ Q ) = ±- J' 2 (l). (75) 
Thus we obtain w(l) = Y (1) J'(l) Since J(l) = 0, 

Y (1) J'(l) = - J(1)Y'(1) + F(l) J'(l) = -W(l) = -1, (76) 
where we have taken the Wronskian from Appendix | A. 2 1 □ 
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5.4 Comparisons of Various Approximate Eigenfunctions 

The outer and inner solutions in our asymptotic expansions of the eigenfunctions may be combined 
into a composite expansion that gives a uniformly accurate approximation in both regions. (See 
Hin91, Section 5.1.8].) This may be formed by adding the inner and outer approximations, 
subtracting the common part in the matching region, and expressing the result as a function of 
the outer variable. This process is trivial at order 0: at y — 1 the inner approximation, and hence 
the common part, vanishes, while at y = 0, the inner approximation is nonzero, but it equals the 
common part, so the uniform approximation is simply the outer solution. It is similarly trivial at 
order 1/3, so let us proceed to order 1/2, where the behavior near y = 1 is interesting. 

Inner and outer solutions are given in Table [TJ and the common part of these expansions in 
the matching region near y = 1 is given by (|71[) or (|72[), 



e 1/2 J'(l)[l- Z\. (77) 

Thus, subtracting off the common part simply cancels the two polynomial terms in the inner 
solution, leaving only the exponential. Therefore, at order 1/2, the composite approximation of 
the eigenfunctions is 

Ui/2(y) = Ay) - e 1/2 X [w(y) - wMe-W*"] . (78) 

Here we have used Lemma [4] to rewrite the coefficient of the exponential to make it obvious 
that ^1/2(0) = 0. Although this boundary condition is satisfied exactly, the derivative boundary 
condition u'(0) = is satisfied only to leading order; specifically W^ 2 (0) = 0(e 1 ^ 2 ). Such loss of 
accuracy in taking derivatives cannot be avoided. 



Fig 3a shows the order-0 outer, the order-1/2 composite, and the order-1/2 inner approxima- 



tions together with the numerical approximation to the first eigenfunction, while Figure 3b graphs 
the derivatives of the order-0 outer approximation and the order-1/2 composite approximation, 
all for e = 10~ 2 . In anthropomorphic terms, Ui/ 2 "attempts" to correct for the fact that the 
zeroth-order outer solution has nonzero derivative at y = 1. Thus, U1/2 is lowered in the interior 
of the interval so that it may approach y = 1 with nearly zero slope. Incidentally, the increase in 
the eigenvalue £ 1 ' 2 Xi/2 is needed to drive U1/2 towards zero more rapidly in the interior. 



In Fig 3a it may be seen that U1/2 differs noticeably from the numerically computed eigen- 
function in the interior of the interval. This difference is corrected by the order-2/3 terms in the 
series, which are driven by the boundary layer at y — 0. (The order-2/3 corrections for clamped 



boundary conditions are the same as for the pinned case, which are calculated in Subsection 6.2 
In other words, for application to the flag problem mentioned in the introduction, it is unwise to 
neglect the order-2/3 corrections to the eigenfunction. 

5.5 Lessons Learned 

Having achieved our initial goal, let us pause to describe patterns in the calculations: Suppose all 
terms u 7 , U 7 , Vy, A 7 of order C(e 7 ) for 7 < a have been calculated, and consider what is needed 
to calculate the terms of order 0{e a ). 



(a) Boundary-layer near y = 0: The equation (261 for U a 



U"" - (XU'Y = RHS, (79) 
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y 0.2 0.4 0.6 0.8 1 

y 

(a) Function. (b) Derivative. 

Figure 3: Comparison of numerical, naive, bulk, and composite approximations of order 1/2 to the first 
eigenfunction, for e = 10 -2 . Both function values and first derivative are shown; in Fig 3a, differences are 
overlaid for clarification. 



has a four-dimensional solution space. However, because J Bi is excluded as unsuited for matching, 
only a three-dimensional space is available for forming U a . This function must satisfy the three 
boundary conditions at X = derived from ( 11|12 1; one expects these boundary conditions to 
determine U a uniquely. 



The particular and homogeneous solutions of (79 1 play different roles. The particular solution 



which has polynomial growth as X — > oo, matches onto derivatives of lower-order terms in the 
bulk series. (For example, the particular solution — XqX in U1/3 matched onto uq^^X).) The 
homogeneous solution, a linear combination A + B J Ai + C^(X), matches onto u a , the bulk 
solution of the same order. Note that as X — > 00 

A + B [ M{x)dx + C^{X) = ClogX + (A + S/3 + C* 00 ) + o(l). (80) 
Jo 

The three-dimensional parameter space is constrained in two directions by matching onto u a . The 
remaining degree of freedom provides the additional flexibility needed to satisfy all the boundary 
conditions. 



(b) Boundary-layer near y = 1: The equation (27 1 for V a , 



V"" - V" = RHS, (81) 

has a three-dimensional space of solutions appropriate for matching. Solutions of the homogeneous 
equation are spanned by 1,Z, e~ z . The constant function matches onto u a , the bulk solution of 
the same order; Z matches onto u Q _i/2(l — e 1 / 2 ^); and e~ z provides the flexibility needed to 



1G 



satisfy all boundary conditions without contributing to the matching. If a > 1, then a particular 
solution of (81 ) will also be needed to match onto derivatives of u a _i, u a _ 3 / 2 > 



The general solution of (81 ) has the form 

V a (Z) = A + BZ + Ce- Z + V P (Z). (82) 

The constants may be determined by satisfying three equations: two come from the boundary 
conditions at Z = and the third arises from matching terms proportional to Z as Z — > oo with 
the derivative of u a _\ji at y = 1. Once V a is determined, matching to u a yields an effective 
boundary condition for u a at y = 1. 

It is interesting to compare matching at the two end points. In both cases, matching constrains 
the inner solution in two directions. At y = 0, both directions relate to u ai the outer solution 
at the order being calculated. In contrast, at y = 1, one direction relates to u a and the other to 

u a-l/2- 



(c) Bulk: In the equation (25 1 for u a , let us split off the term on the right proportional to the 



yet-to-be-calculated coefficient X a : 

(yu'Y + X u = —X a Uo + R. 
The general solution of this equation has the form 

u a — —\ a w + u p + aJ + bY, (83) 



where w is the solution of (581 and a,b are arbitrary. If it weren't for the complications arising 
from the logarithmic behavior of Y near y = 0, the remaining steps would be extremely simple: 
Matching u a to U a at y = provides two equations, which we may use to determine a and b in 



(83 1. Then the effective boundary condition from matching at y = 1 provides a linear equation 
for A Q , completing the calculation to order 0(e a ). These simple ideas suffice until we encounter 
order a = 1; even then, the complications are rather mild. 

(d) Matching: To match at y = we need 

^ [u,{e^X) - U 7 (X)] = o(s a ) (84) 

for an appropriate range of X. Since all terms of order less than a have already been matched, 
we may focus only on terms of order exactly a. In the inner series, only U a contributes a term of 
order exactly a: 

Y.^U^X)^ ... + e a U a (X) + o(e a ) 

where . . . indicates terms of order less than a, which are not relevant for calculating the 0(e a )- 
solution. By contrast, in the outer series, in addition to u a , derivatives of lower-order terms in the 
series also contribute terms of exactly this order. Of course these lower-order terms in the outer 
solution have already been determined; typically they are matched by the particular-solution part 
of U a , the arbitrary constants in the homogeneous solution playing no role. The "bleeding" of 
lower-order terms into the order-a matching is also responsible for the fact that the matching 
interval shrinks as a increases. 

Similar considerations apply to matching at y = 1. 
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We have everywhere performed matching in terms of the inner variable. It is possible to use 
the outer variable instead, but in our opinion the calculations are less clear: Specifically, the terms 
needing to be matched at order a are precisely those that, when expressed in terms of the inner 
variable, are proportional to e a . 



6 Asymptotics for the Pinned Case 
6.1 The Low-Order Solution 

We now consider pinned boundary conditions at y — 1 with the same conditions at y = 0: 

u"(0)=Q, u"'(0)=0, u(l)=0, u"(l)=0, (85) 

plus the normalization u(0) = 1. We calculate the terms in the asymptotic series through the 
first nontrivial correction to the eigenvalue — in this case order one. The results are summarized 
in Table [2 





1/3 
1/2 
2/3 
5/6 
1 



Ag 

A 







J(y) 






6Ai'(0) 





J(y) 



^x 2 



u* 

1 

-X X 


2Ai'°(o) Jo Aifc) dx 




Vg 




-J'(1)Z 





(1111 



(1121 



(94) 



(97) 



Table 2: The coefficients of the asymptotic expansion for the pinned case. The last row contains the 
equation numbers where the terms, too long and cumbersome for the table, can be found. 



Exercise 2. Derive the first 3 rows of Table^ 

These first few orders are very similar to the clamped case. 

6.2 The e 2 / 3 Correction Term 

(a) Boundary-layer near y — 0: Using the values of and from Table [2j the reduced 



Therefore 



equation (26) for [/ 2 /3 has the inhomogeneous term — AqX, with particular solution U p — -^-X z 



the relevant family of solutions is 



= ^x 2 



x 



A + B Ai(r)dr + C ■ 
Jo 



The boundary conditions C/ 2/3 (0) = f7£' /3 (0) = U^ 3 {0) 



determine that A 



(86) 

C = and 



B = — 



2Ai'(0) 



So we have verified the table entry for U 2 /3- 
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(b) Boundary-layer near y = 1; Equation (27) together with the two boundary conditions at 
y = 1 imply that V2/3 has the form 



V 2/3 (Z)=A(l-Z-e- z ). 



(87) 



(c) Bulk: Because vanishes, the ODE (251 for w 2 /3 has only the inhomogeneous term 
-A2/3U0 The general solution of this equation is 

"2/3(2/) = -A2/3WG/) + aJ(y) + bY(y), 



where w(y) is the function that solves (58 1. 



(d) Matching: At y = we need equation (84 1, with a — 2/3, to hold for X in the range 
e~ p < X < £-« where < p < q < 1/3. We take p = 0,q = 1/9 so that (e 1 / 3 ^:) 3 = o(e 2 / 3 ), and 
we focus only on terms of order exactly 2/3, using ellipsis to represent terms of lower order. From 
(c) and using the fact that w(0) = 0, 



£~ t u 7 (e 1/3 X) = 



-2/3 



7 <2/3 



<(0) 



X 2 + a + blog{e 1/3 X) 



+ o(£ 2/3 ) 



On the other hand, reading U2/3 from Table[2] and using (fl20|) from the Appendix |A. 1 



7 <2/3 



-2/3 



A v 2 



x z - 



A 2 



6Ai'(0) 



(89) 



We see from the Taylor series (126) that the quadratic terms in (88) and (89) agree, as expected. 
Matching the new terms at order 2/3 gives 



\2 



6Ai'(0) 



6 = 0. 



(90) 



At y = 1 we match for Z in the range e p -C Z -C £ q where we take p — 0,q — 1/6 so that 
(e 1 ^ 2 Z) 2 — o(£ 2 / 3 ). Again we focus only on terms of order exactly 2/3. Now 



E £ 1 u 1 {l-e 1/2 Z) = 

7 <2/3 



+ £ 2 / 3 u 2/3 (l) + (£^ 3 ); 



2/3N 



no derivative of the bulk solution bleeds into the £ 2 / 3 term since there is no term of order 1/6. At 
the same time, by ((87|) 



£ ^u,{z) 

7 <2/3 



+ Ae 2 / 3 (l-Z)+o(e 2 ^). 



Thus matching requires A — and 1*2/3(1) = 0; the latter may be written out as 

-A 2/3W (l)+ a J(l) = 



where a is given by (90). Of course J(l) = and w(l) ^ 0, thus it follows that A 2 / 3 = and m 2 / 3 
is as given in Table 27^ 
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(e) The composite approximation: As discussed in Subsection 5.4 the composite approximation 
is given by the sum of the outer and inner approximations minus the matching terms. Taking the 
outer and inner solutions from Table [2] and the matching terms near y 



2/3 



\2 



6Ai'(0) 



J(y) - 3 



from (|89|, we find 
AIM* 



(91) 



Here the integral term results from a convenient cancellation: in the matching layer, as X —¥ oo, 
the integral over (0, X) tends to the integral over (0, oo), so the inner solution minus the matching 
terms simplifies to the integral over [X, oo). 



Figure 4a compares the numerical solution and two other approximations to the third eigen- 
function: the naive approximation and the order-2/3 composite approximation. The graphs of 
all three functions resemble that of Figure 1, and to visual accuracy they coincide; therefore we 
have plotted only differences. Figure [4b] graphs the second derivatives of these three approxi- 
mate cigenfunctions. In all plots e = 10~ 5 ; because we are examining the third eigcnfunction, we 
need a smaller e than in Subsection 5.4 in order for the asymptotics to be meaningful — see Sec- 



tion [7] The order-2/3 composite approximation captures most of the boundary-layer behavior of 



the eigenfunction. In particular, observing in the Taylor series (126) that J"{0) = we compute 
that 

U" ,(0) = ^-£ 2/3 ; 

2/3 v ; 12Ai'(0) 

thus, the second derivative vanishes to leading order. 




(a) Differences. (b) Derivatives. 

Figure 4: Comparison of (a) function values and (b) second derivatives of order-2/3 approximations to 
the third eigenvalue, for e = 1CP 5 . Since the function values are too close to discern visually, various 
differences between functions are shown instead. 
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6.3 The e 5/6 Correction 

Exercise 3. Show that all the terms in this order vanish, as shown in Table^ 

6.4 The e 1 Correction 

This order of the calculation is pedagogically interesting because a logarithmic term enters the 
series. 



(a) The y = boundary: Equation (26 1 for Ui is 



U'i" - (XU[)> = |x 2 - [ X Ai(x)dx. (92) 



o 

This looks quite bad due to the integral of Ai in the RHS, but in Exercise [4] you are guided through 
a proof that 

^--£* + 5^jf ( *-' )JU(,>fc < 93 > 



is a particular solution to ( 92 ) , and in Exercise [5] you are asked to apply boundary conditions to 
conclude that 

U,(X) = G(X) - | £ AKx)dx + (94 ) 



Exercise 4. Find a particular solution to (|92j : 

1. Show that the derivative of J* (X — x)Ai(x)dx equals Ai{x)dx. 

2. Use this when applying the LHS of ODE ( |92[ ) to J (X— x)Ai(x)dx and obtain — J Q Ai(x)dx. 



3. Conclude that G(X) in (93) is a particular solution of ( |92[ ) that grows only polynomially at 
large values of X. 

Exercise 5. Find the y = boundary layer solution U\: 

1. Using Airy's differential equation, show that 

G(0) = 0, G"(0) = |^, G"'(0) = |. (95) 

2. The family of possible solutions for U% is 

cX 



U 1 (X) = G(X) + A + B I ki{x)dx + C^{X). 
Jo 



(96) 



Using what you just found about G, show that (94) satisfies the boundary conditions E/i(0) 
U['{0) = U["{0) = 0. 
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(b) The y = 1 boundary: The boundary layer equation is V{"' — V{' = J'(l) whose family of 
possible solutions is 



Vi = - 



J'(l) 



A + BZ + Ce 



-z 



The boundary conditions V(0) = V"(0) = determine A and C so that 

1 



V 1 = -J'(l) 



-Z 2 



1 



BZ, 



(97) 



and by matching to the bulk solution (below), it is found that B = 0. 



(c) Bulk: The bulk equation for u\ is familiar, albeit with a new term on the RHS: 

(y<)' + A oUl = -Aiuo + (98) 

Its general solution is 



= -Aiw(y) + + aJ(y) + bY(y), 



where w(y) is defined by (58 1 and v(y) is the solution to 

(yvj + X v = J""(y), v(0) = 0. 
(Existence and uniqueness of v(y) provided by Lemma [T]) 



(99) 
(100) 



(d) Matching: Near y = we need (84) to hold for a range of X such that e~ p <C X <C e~ q 
where < p < q < |. We take p = 0, g = 1/12 so that (e 1 / 3 ^) 4 = o(e 1 ). Using an ellipsis for 
terms of order less than 1, we have for u 



^e^ Ul {e 1/z X) 

7<1 



^) X s_>MMx + a. 



6 



6Ai'(0) 



1 



log e + log X 



+ o(s 1 y, (ioi) 



the cubic and linear terms come from derivatives of uq and U2/3 respectively. At large X, the 
behavior of U\ is 



7<1 



where 



C' 



Ai(0) 



+ o(e 1 ) 



0.06855. 



6(Ai'(0)) 2 ' 3 

Here the 1/2 in comes from a term in G(X), using the relation 

xAi(x)dx = -Ai'(0) 



(102) 



(103) 



obtained from the Airy equation. It is readily seen from ( 126 ) that the cubic and linear terms 
of (101 1 and (102 1 match, and the logarithmic terms will match if b = Aq/3. However, it is not 
possible to match ( 101 1 and ( 102 ) with a coefficient a that is independent of e. To fix this problem 
we propose to augment the u-series with a logarithmic term: i.e., to replace eu\(y) by 

e[(bgs) ui{y) + ui{y)\. (104) 

In Exercise [6] below we ask the reader to show that: 
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No such log terms are possible in the asymptotic series for either boundary layer or for A. 



• The equation for u± is just the homogeneous version of (25 1 so that for some coefficients d, b 

Ui{y) = aJ{y) + bY{y). (105) 



Therefore in matching, the RHS of ( 101 ) should be replaced by 
>'(0) X 3 AgJ'(O) 



+ e 



(i 



6Ai'(0) 



X + (a + aloge) + (b + bloge) - logs + log A 



1 



Matching between ( 102 ) and ( 106 1 is now possible if and only if 

a = C ao , b = \ 3 /3, a = -6/3 = -A 3 /9, 



6 = 0. 



+ o(e 1 ). (106) 



(107) 



Near the boundary at y = 1 we match for Z in the range 1 <C Z <^ e x / 6 so that both 
(e 1/2 Z) 3 = o(e 1 ) (to allow neglect of u' ") and e 2/3 ■ e 1/2 Z = o(e x ) (to allow neglect of u' 2/3 ). As 
Z— >oo, the u — series has the asymptotic behavior 



+ £ 



ug(l)^- + Ul (l) 



o(e); 



(108) 



mercifully the contribution of u\, which is proportional to </(l), vanishes. From (97) the 1^-series 
has asymptotic behavior 



-J'(l) 



2 



1 



BZ 



o(e). 



(109) 



According to Exercise [7] below, the quadratic terms match; for the linear terms to match we need 
B = 0; and for the constant terms to match we need 



Ul (l) = - J'(l). 



Thus manipulating ( 99 ) , we deduce that 
Ai 



„(l) + £y(l) + J'(l) J Q L J""(r)J(r)d7 
= Aq- 



1 A* 



w(l) 



3 [J>(1)Y 



Ao, 



(110) 



(111) 



the latter equality following from the proof of Lemma [T] from Lemma |4j and ( 76 ) . We may 
regard the three terms in (111) as contributions from the fourth derivative in ([8]), the boundary 
conditions at y = 0, and the boundary conditions at y = 1, respectively. (Numerically, we find 
that Ai w 4.4280, 1887.2 and 44403 for the first three eigenvalues.) Finally, we have 



A 3 A 3 
(k>ge)fii(i/) + u x (y) = -Ai«;(y) + w(y) + (C^ - loge) J(y) + f (y) 



(112) 



where is given by (103). 

As a check on our calculations, in Figure [5a] we present a log-log plot of the error in the two- 
term approximation Ao + Ai£ to the eigenvalue. The remaining error in the approximation seems 
to have a slope of 4/3. 
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(b) Clamped case. 



Figure 5: A log- log plot of the error of the two-term approximation to the eigenvalues, together with a 
line of appropriate slope for visual reference. 



Exercise 6. Suppose that all of the series — for u, U , V , and A — contain terms of order eloge, 
with the relevant functions identified by a "hat. " 

1. Deduce that U\{X) = because it satisfies a homogeneous differential equation of the type 
(26 1 with homogeneous boundary conditions. 



2. Deduce from the (homogeneous) equation for V\, of the type (27), and the (homogeneous) 
boundary conditions that V\(Z) = BZ for some constant B. Show by matching with the 
lower-order parts of the inner solution that B = 0. 



3. Observe that U\{y) will satisfy an ODE of the type (25) 

(yu'iiv))' + ^oUi{y) = Aiu - 



Deduce from this equation plus matching that Ax = and U\ has the form (105). 



Exercise 7. Using the fact that J(l) = 0, deduce from ODE (25) for uq that 

J'(l) = - J"(l). 

We invite the ambitious reader to calculate the additional terms in asymptotic series for the 
clamped case through 0{e): The order-2/3 terms are identical and the order-5/6 terms have 
one minor difference, while there are some significant differences at order-1. For example, the 
correction for the eigenvalue is given by 



Ai = 



v(l) + Agr(l) + %Y (1) - f J'(l) + X Q w'(l) 
w(l) 



(113) 
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where r is defined by 

(yr'Y + X r = w, r(Q) = 0. 
Numerically, Ai sa 6.4581, 1900 and 44435 for the first three eigenvalues. 



7 Closing Remarks 

As mentioned in the introduction, the aim of this paper was to understand the effect on the 
eigenf unctions for ^ of a small, but nonzero, bending resistance. We have found that, for small 
epsilon, the order-2/3 composite approximation (91 1 provides a reasonable correction to the naive 
eigenfunctions. However, the phrase "small epsilon" needs clarification. 

Our notation so far hides the fact that problem pj has an infinite sequence of eigenvalues 
\( n \e) that tend to infinity as n — > oo. For any fixed n, \( n '(e) has an asymptotic series in e with 



(n) - - 

leading order A characterized by (17). We ask, for a given n, how small must e be to get into 
the asymptotic range. A minimal requirement is that the order-2/3 correction to the eigenfunction 
be small compared to unity: i.e., 

e 2/3 X 2 < 1. (114) 



To make this more quantitative, we invoke the asymptotic approximation (122) of the Bessel 
function to deduce that, asn->oo 



An) 



1n2 I" 
A> 4 



Thus, (114) requires that 



e < n u . 

In other words, the meaning of "small epsilon" depends heavily on n. 



(115) 



(116) 



In fact, for some purposes, even ( 116 ) is not sufficiently restrictive. Let us define the asymptotic 
region by requiring that, in Figure \A say in the clamped case, the graph of A("'(e) - A ™ } has 
converged to a line of slope 1/2; in symbols, 



|A (n) (e)-A 



(n) 



.1/2 \(n) 
■ A l/2' 



(117) 



We may estimate the LHS of (117) by the next term in the asymptotic series, eA$™\ and it may 



be shown by estimating the terms of (113) that X^ = 0{n 7 ). Since A^ = A o"^ anc ^ by (115) the 



>(«) 



v(") 



latter is 0(n 2 ), formula (117) is equivalent to en 7 <C e 



r l/2_2 



or e <C n 



1/2 
-10 



The restrictions in applying small-e asymptotics are so severe that often problem ([8j) is more 
accurately approximated by treating e as a large parameter: i.e., regarding p| as a perturbation 
of a fourth-order equation by a second-order operator. This point of view may be developed 
systematically to explain the straight-line behavior in the range 10 -2 < e < 1 in Figure [2] and 
predict that this behavior will continue indefinitely as e increases. 
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A Properties of Special Functions 

The properties of the Bessel and Airy functions in this Appendix are taken from Chapters 9 and 10 
of |AS64] . respectively. 



A.l The Airy Functions, Ai and Bi 

The Airy functions, Ai and Bi, are the standard solutions to the Airy equation u" — xu = 0. The 
Wronskian of these solutions is constant: 

M(x)Bi(x) - Ai'(a;)Bi(a;) = -. (118) 

7T 

The Airy functions have exponential behavior for large positive x: 

_2j.3/a l x 3/2 

Ai(a)~ 6 ' Bi(x)~-^™. (119) 

The definite integral of Ai converges and is known explicitly |AS641 10.4.82]: 

1 

Ai(a:) = -. (120) 



Incidentally, 



Ai/ (°) = -3W(l73)- (121) 



A. 2 Bessel Functions of Order Zero 

The two standard solutions to the zeroth-order Bessel Equation u" + ~u' + u = are denoted Jo 
and Yq. Despite the singularity of Bessel's equation at x = 0, Jq(x) is differentiable near zero and 
in fact its power series converges for all x. By contrast, as x — > 0, 

Y (x) = 2 (log(ix)+ 7 ) +0(x 2 \ogx) 
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where 7 = 0.577215665 ... is Euler's constant |AS64| 9.1.13]. For large x, these functions have the 
asymptotic approximation 



Jo(x) ~ \ — cos(x - 7T/4), 

V irx 

Y o(x) ~ \ — sin(x-7r/4) 

V TTX 

In particular, Jq has an infinite sequence of positive real zeros, beginning with 



(122) 
(123) 



2.40482, 5.52007, 8.65372, .... (124) 

In the paper we need linear combinations of Bessel functions with a slightly different argument 
as defined in ( 29 1 and ( 30 ) . The Wronskian of J and Y may be computed from the fact that the 

1/ttx and making the substituti 

W = J(y)Y'(y)-Y(y)J'(y) = 



Wronskian of Jo and Yq equals 1/ttx and making the substitution x = 

1 

y' 

J has the power series representation 

J{y) = 1- X y + 



2! 



( 

{ 3! 



(125) 



(126) 



This may be obtained by substitution into the power series for Jq(x) or by computing coefficients 
recursively from the ODE (14). 



B Proof: The Boundary-Layer Function \I/ 



Proof of Lemma^ (Existence) Equation ( |41[ ) has an irregular singular point at X — 00 but is 
regular everywhere else. Near the singular point, this equation admits a formal series solution 



V(X) ~ \ogX 



OO 

J^c k X~ sk 

k=l 



(127) 



where C\ — 2/3 and subsequent coefficients may be determined recursively. Therefore there is a 
solution of (41), an entire function of X, that has the asymptotic series (127) as X — > 00, and in 



particular, the asymptotic behavior in the lemma. By subtracting off a multiples of U^-\ we 
can also satisfy the boundary conditions at X — without affecting the logarithmic behavior as 
X^oo. 



(Uniqueness) Now suppose 'J is another solution of (41) that also satisfies the conditions of 
the lemma, and let $ = * - Then <I> solves ([36]) and W{X) = 0(1) as X^oo. $ must be a 
linear combination of {U^\ i = 1,2,3,4}, where these functions are defined in Subsection 4.4 
However, and both grow more rapidly than 0(1), and so does any nontrivial linear 

combination of them. We conclude that <f> = cif/t 1 ) + C2U^ for some constants ci,C2, and since 
$(0) = $"(0) = 0, we find that c x = c 2 = 0. □ 

Figure [6] shows a numerically computed graph of ty(X). It was found by applying a boundary- 
value solver (Matlab's bvp5c) to (41), specifying two boundary conditions at X — and an 



estimate for the first derivative at X — 20 obtained from the first 6 terms in (127) 
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